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ABSTRACT 

We present optical spectroscopic measurements of the eclipsing High Mass X-ray Binary 
XMMU J013236. 7+303228 in M33. Based on spectra taken at multiple epochs of the 1.73 d bi- 
nary orbital period we determine physical as well as orbital parameters for the donor star. We find 
the donor to be a B1.5IV sub-giant with effective temperature is T = 22, 000 — 23, 000 K. From the 
luminosity, temperature and known distance to M33 we derive a radius of i? = 8.9 ± 0.5 i?0. From 
the radial-velocity measurements, we determine a velocity semi-amplitude of Kopt = 63 ± 12kms~^. 
Using the physical properties of the B-star determined from the optical spectrum, we estimate the 
star's mass to be Mopt = 11 + 1 Mq. Based on the X-ray spectrum, the compact companion is likely 
a neutron star, although no pulsations have yet been detected. Using the spectroscopically derived 
B-star mass we find the neutron star companion mass to be Mx = 2.0 + 0.4 M©, consistent with the 
neutron star mass in the HMXB Vela X-1, but heavier than the canonical value of 1.4 Mq found 
for many millisecond pulsars. We attempt to use as an additional constraint that the B star radius 
inferred from temperature, flux, and distance, should equate the Roche radius, since the system ac- 
cretes by Roche lobe overflow. This leads to substantially larger masses, but from trying to apply 
the technique to known systems, we find that the masses are consistently overestimated. Attempting 
to account for that in our uncertainties, we derive Mx = 2.2^Q gMo and Mopt = 13 + 4Mo. We 
conclude that precise constraints require detailed modeling of the shape of the Roche surface. 
Subject headings: none 



1. INTRODUCTION 

The range of possible neutron star masses depends 
on many factors, such as the initial mass of the pro- 
genitor's stellar core, the details of the explosion (in 
particular mass accretion as the explosion develops), 
subsequent mass accretion from a binary companion, 
and the pressure-density relation, or equation of state 
(EOS), of the neutron star matter. On the low- mass 
end, producing a neutron star requires the progenitor's 
core to exceed the Chandrasekhar mass, which depends 
on the uncertain electron fraction. Theoretical models 
place this minimum mass in the range Mcorc ^ 0.9 - 
1.3 Mq (jrimmcs ct al. 1996;). The largest possible neu- 
tron star mass depends on the unknown physics deter- 
mining the EOS - for example whether kaon condensates 
or s trange matter can from in t he interior (See, for exam- 
ple, |LattimerXPra^i2|200l). The highest -mass neu- 
tron star to-date with an accurate measurem ent weighs 
in at 1.97 + 0.04 Mq ("De morest et al.ll2010D . which al- 
ready rules out the presence of e xotic hadronic matter 
at the nuclear satura tion density ([Demorest et al.|[2010l : 
iLattimer eralll2010( ). 

Testing the predictions of supernova models, binary 
evolution models, and finding objects at the extremes 
of the mass spectrum require determining neutron star 
masses in a variety of systems with differing progenitor 
masses and evolutionary history. Neutron stars accom- 
panying either a high-mass star or another neutron star 
are thought to have accreted little to no matter over their 

^ Cahill Center for Astrophysics, California Institute of Tech- 
nology, Pasadena, CA 91125, USA 
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lifetimes. In contrast, neutron stars in low-mass X-ray 
binaries and millisecond pulsars, typically in close orbits 
around a white dwarf, have undergone extended accre- 
tion periods that will make the current mass exceed that 
at birth. Different types of binaries will also have differ- 
ent average neutron star progenitor masses. 

High-mass X-ray Binaries (HMXBs) - binaries contain- 
ing a neutron star and a massive (~ 20 Mq) companion 
- are particularly interesting systems in which to pursue 
mass measurements. In most cases the neutron star pro- 
genitor will have been more massive than the observed 
donor star, yielding a relatively high-mass pre-supernova 
core. Furthermore, the NS mass will be close to the birth 
mass, since even for Eddington rates < O.IMq can be 
accreted in the ~ 10^ yr lifetime of the OB companion. 
Indeed, among the five HMXB with reaso nably secure 
masses, one (Vela X-1) has M = 1.8 AfQ ()Barziv et al.l 
l2001t[Quaintrell et al.ll2003l ). indicating that this neutron 
star may have been born heavy. 

Determining NS masses in HMXBs is, however, diffi- 
cult. In compact object binaries (NS - NS, NS - white 
dwarf), highly precise mass measurements can be ob- 
tained from relativis tic effects like the precession of peri- 
astron (iFreire et aLll2008i) or measurement of the Shapiro 
delay (jPemorest et al.ll2010D . In HXMBs, however, accu- 
rate mass measurements are limited to eclipsing systems 
where orbital parameters for both the NS and its stel- 
lar companion can be measured. For the NS this is done 
through X-ray or radio pulse timing, and for the compan- 
ion through radial-velocity measurements derived from 
doppler shifts in the stellar lines. In the event pulsations 
are not detected, the NS mass can still be determined if 
good spectra are available to estimate the mass of the 
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optical component. In rare cases where the distance to 
the binary is known, this provides an independent con- 
straint on the physical scale of the system - for exam- 
ple by calculating the absolute magnitude of the compo- 
nents. However, calculating masses from such constraints 
is model-dependent. 

In this paper we present optical spectroscopic 
measurements of the donor star in the eclips- 
ing HMXB XMMUJ013236. 7-^303228 using the 
Low Resolution Imaging Spectro grap h on t he 10 m 
Keck-I telescope (LRIS; Oke e t al.l 119951 ) aimed 
at determining the mass of the compact compan- 
ion. XMMU J013 236.7-I-303228 was discovered by 
iPietsch etall (|2004D in their XM M-Newton survey of 
M 33 . In follow-up observations, iPietsch et al.l ()2006l 
120091) identified it as an eclipsing High Mass X-ray 
Binary with a 1.73 d period. The X-ray spectrum is 
hard, and the shape implies that the compact object is 
a neutron star. However no pulsations were detected 
in the X-ray data, s o a black hole cannot b e ruled out 
(IPietsch et all I2OO90 . iShporer et all (|2006D discovered 
an optical counterpart (Figure [T]) which shows vari- 
ability consistent with ellipsoidal modulation of the 
OB star. Given high quality spectra we are able to 
obtain a spectroscopic mass for the donor, and therefore 
determine the compact object mass. Using the known 
distance to M33 combined with the B stars luminosity 
and temperature we derive a physical radius, which we 
equate with the Roche radius based on the observation 
that accretion is occurring via Roche lobe overflow. This 
provides an additional orbital constraint that we use to 
independently estimate the compact object mass. 
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Fig. 1. — A V-band CFHT image showing the optical counter- 
part to XMMU J013236.7+303228, located at a = 0l'"32™36.=94, 
5 = -|-30°32'28."4 (J2000). The image was obtained from 
CFHT online data archiv e, and a WCS was added using astrome- 
try.net l ILang et al.ll20ia) . 



We observed XMMU J013236. 7-1-303228 on UT 2009 
October 16 and October 17, with the Low Res- 
olution Imaging S pectrograph on the 10 m Keck- 
I tel escope (LRIS; IQke et all 119951), with upgraded 
blue (IMcCarthv et al.lll998l: ISteidel et al.ll2004[) and red 
cameras (jRockosi et al. I I2Q10I) . covering a wavelength 
range from 3,200 A -9,200 A. We set up LRIS with the 
600/4000 grism on the blue side and the 600/7500 grat- 
ing on the red side, to get dispersions of 0.6A/pix and 
0.8A/pix respectively (Tabled]). To maximize stability 
of the spectra, we used the "stationary rotator mode" , 
where the instrument rotator was held fixed near zero 
degrees rather than tracking the parallactic angle while 
observing. Atmospheric dispersion was compensated for 
by the ADC (Atmospheric Dispersion Corrector). We 
acquired a total of 28 spectra of the target, with expo- 
sure times ranging from 300-1800 seconds. The spec- 
trophotometric standard EG 247 was observed for flux 
calibration. 

The data were reduced in IRAF0- The spectra were 
trimmed and bias subtracted using overscan regions. No 
flatfielding was applied. Cosmic ra ys were rejected us- 
ing L. A. Cosmic (|van Dokkum|[200ll ). At mospheric lines 
are s table to tens of meters per second (jFigueira et al.l 
|2010( ). so the wavelength solution for the red side was 
derived using sky lines for each image. For the blue side, 
the wavelength solution was derived from arcs taken at 
the start of the night. The spectra were then rectified 
and transformed to make the sky lines perpendicular to 
the trace, to ensure proper sky subtraction. The wave- 
length solutions for arcs taken at various points during 
the night are consistent with each other to a tenth of a 
pixel, with only an offset between different arcs. We cor- 
rected for this offset after extracting the spectra, by using 
the 5577.34 A [O I] line. The spectra were extracted with 
APALL, and flux calibrated with data for EG 247 and the 
standard IRAF lookup tables. We further tweaked the 
flux calibration by using a EG 247 model spectrum from 
the HST Calibration Database Archivfl We used just 
one standard spectrum per night, and enabled airmass 
correction in IRAF during flux calibration. 

The flnal spectrum is shown in Figures [5] & El The 
signal-to-noise ratio per pixel is > 10 for most blue side 
spectra (Table [Ij. The first night, we experienced some 
tracking issues with the telescope, so the target did not 
remain well centered on the slit at all times. A similar 
problem was experienced for the first exposure on the 
next night, where the object was at a high airmass. In all 
following discussions, we reject some such spectra based 
on observing logs, and spectra with signal-to-noise ratio 
per pixel < 10 (Tabled]). 

3. DONOR STAR PARAMETERS AND ORBIT 

We determine the best-fit stellar parameters and or- 
bital solution using an iterative technique. First, we es- 
timate a spectral type for the primary (donor) from in- 
dividual spectra. We use appropriate spectral templates 
to calculate the orbital solution (Section [3. 2p . Next, we 
shift spectra to the rest frame and combine them to get 
a higher quality spectrum. We calculate stellar param- 
eters from this combined spectrum, and use a template 



2. OBSERVATIONS AND DATA REDUCTION 



^ http://iraf.noao.edu/ 

* ftp:/ /ftp. stsci.edu/cdbs/current_calspec/ 



Constraints on the Compact Object Mass in the Eclipsing HMXB XMMU J013236.7+303228 in M 33 



3 



TABLE 1 

Details of individual exposures 



Image Name 


MJD 


Seeing'' 
arcsec 


SNR 


— 

Heliocentric Velocity 
km 







Notes= 


buyiuiD-UUDz 


5512U.27D 


1.0 


5.2 


-106 it 24 


0.96 


Kejected: 


Logs, Foor bJNH 


d09101d_00d3 


55120. 28o 


1.0 


21.7 


-49 it 5 


1.05 


Selected 




v,nni ni a. c\c\ci.A 
bUyiUiD_UUD4 


r: cr 1 on Qnc; 

ooizU.oUo 


1.0 


23.0 


-32 it 5 


1.08 


Selected 




buyiuiD_uuDD 


c cr 1 on 1 oo 


1.0 


22.8 


-57 it 5 


1.20 


Selected 




v /"in 1 ni c c\c\ o n 

buyiuiD-Uuou 


c cr 1 on Qn /I 


1.0 


22.7 


-101 it 6 


1.36 


Rejected: 


Logs 


bUyiUlb.UUol 


c c 1 on ,110 

55121). 418 


1.0 


8.7 


-109 it 15 


0.99 


"r> ^ \ ^ -J 

Kejected: 


Foor bNH 


buyiuiD-Uuoz 


c cr 1 on -IOC 

5512U.425 


1.1 


20.1 


-50 it 6 


1.16 


Selected 




bUyiUlb-UUoY 


c cr 1 on A 

5512U.457 


1.1 


13.7 


-47 it 9 


1.09 


Selected 




"Lnm ni n nnoo 


c cr 1 on A'~? A 

5512U.474 


1.4 


4.2 


-65 it 30 


1.14 


"r> ^ J 

Kejected: 


Logs, Foor bJNH 


buyiuib-Uuyo 


c c; 1 on c; m 
0512U.5U1 


1.2 


18.3 


-103 it 6 


1.08 


Rejected: 


Logs 


buyiUib-Uuyy 


cc c; 1 on c; o /i 
00l2U.024 


1.1 


19.0 


-68 it 6 


1.02 


Selected 




bUyiUib.UlUb 


c cr 1 on cr rT\ 

5512U.5d(J 


1.2 


15.5 


-131 it 8 


1.22 


Selected 




v,nni m a. ni no 
buyiUib_UiUo 


c cr 1 on cr o^; 

o512U.5ob 


1.1 


11.4 


-119 it 11 


1.00 


Selected 




v /"in 1 ni c m nn 

buyiuiD-Uiuy 


c cr 1 on cr n^v 

5512U.5y7 


1.1 


7.6 


-162 it 19 


1.09 


"r> ^ \ ^ -J 

Rejected: 


Foor biNH 


bUyiUlb-UiiU 


c: e; 1 on i^nQ 
5512U.bUo 


1.2 


6.3 


-izf5 ih 24 


1.03 


Rejected: 


Logs, Foor oiNrt 


v,nni m i nnc^rr 


c; cr 1 o 1 ovn 
55121.2 (M 


0.9 


11.3 


-42 it 11 


0.99 


Rejected: 


Logs 


hDQi ni 7 on^Q 

\J\JU 1-\J 1- 1 —\J\J\JZj 




0.7 


18.8 


-96 it 6 


0.94 






b0yi017.0060 


55121.327 


0.9 


iy.6 


-79 it 6 


0.93 


Selected 




b0yi017_0068 


55121.367 


0.8 


22.5 


-75 it 5 


0.93 


Selected 




b0yi017.0069 


55121.389 


0.8 


22.5 


-61 it 5 


0.97 


Selected 




b0yi017.0071 


55121.419 


1.0 


iy.5 


-77 it 6 


0.99 


Selected 




b0yi017.0078 


55121.457 


1.0 


iy.7 


-57 it 6 


0.95 


Selected 




b0yi017.0079 


55121.478 


1.1 


16.9 


-27 it 7 


1.03 


Selected 




b0yi017.0081 


55121.507 


1.1 


18.0 


-46 it 7 


0.98 


Selected 




b0yi017.0087 


55121.541 


1.1 


17.2 


-44 it 7 


1.02 


Selected 




b0yi017_0089 


55121.567 


1.1 


15.1 


-53 it 8 


1.03 


Selected 




boyioi7.ooyo 


55121.589 


1.1 


8.7 


-70 it 13 


0.89 


Rejected: 


Poor SNR 


boyioi7.ooyi 


55121.600 


1.1 


8.4 


-32 it 15 


0.88 


Rejected: 


Poor SNR 



Note. — LRIS was set up with a \" slit, the D560 dichroic and clear filters on both red and blue arms. To 
maximize stability, we used the stationary rotator mode. The 600/4000 grism gives a dispersion of 0.61 A/pix on 
the blue side. On the red side, we configured the 600/7500 grating at 27°. 70 (central wavelength 7151A), and 
get dispersion of 0.80A/pix. 

^ Seeing was measured as the median FWHM of Gaussians fitted to the trace of the blue side spectrum at several 
points. 

Errors quoted here do not include an additional 15kms~^ error to be added in quadrature, due to motion of 
star on the slit. 

Spectra were not included in the final analysis if either the observing logs mentioned that the star had moved 
from the slit, or if the signal-to-noise ratio (SNR) per pixel in the extracted spectrum was below 10. 



spectrum with these refined parameters to recalculate the 
velocities. In this section, we describe the final iterations 
of both these steps. 

3.1. Stellar parameters 

Based on p hotometry of XMMU J013236. 7+303228, 
iPietsch et alT (j2009[) estimate that the companion is a 
10.9 Mq object with T^s = 33000 K and \og{g) = 4.5, 
with = 2.4 for their best-fit model. They then as- 
sume a distance of 795 kpc to M33, and calculate that 
the star has an absolute magnitude My ~ —4.1 and the 
line of sight extinction is Ay — 0.6, so derive a stellar 
radius of 8.0 i?©. 

We deduce the spectral type by comparing our spectra 
to iWalborn fc Fitzpatrickl ()1990t ) and the Gray spectral 
atlafl The absence of He II lines (Figure [2]) implies a 
spectral type later than O, while the relative strengths 
of the Mg II 4482 A/He I 4471 A lines point to a spec- 
tral type earlier than B3. The strength of He I lines, 
and a weak feature at 4420 A indicate a spectral type 
around Bl, for main sequence stars. A bump blueward 
of H8 3889 A is characteristic of spectral type B2. The 
weakness of C III 4650 A, and the relative strengths of 

^ http:/ /ned. ipac.caltech.edu/level5/Gray/frames. html 



C III/O II near 4650 A refine the spectral type to be- 
tween Bl and B2 for both dwarfs and giants. Finally, 
the weakness of Si IV RS 4101 A, gives a spectral type of 
B1.5. To determine the luminosity class, we note that the 
O II 4415-4417 A and Si III 4552 A lines are present, 
but are weak as compared to He I 4387. We conclude 
that the donor is a B1.5IV star, with rough uncertainties 
of 0.5 spectral subclasses and one luminosity class. 

Tabulated values of stellar parameters are usually pro- 
vided for luminosity class HI and V stars. For B1.5V 
stars, T « 23000 K, log(5) = 4.14, and My = -2.8. 
For B1.5III s tars T » 22000 K, log(g) = 3.63, and 
My = -3.4 (|Co3[2000f) . The B1.5IV target wiU have 
values intermediate to these. This inferred tempera- 
tu re is significan t ly co oler than T 33000 K reported 
bv IPietsch et al.l (|2009f ). However, the absence of He II 
lines at, e.g., 4541 and 4686 A, clearly exclude such a 
high temperature. 

The colors of a star depend on its temperature and 
surface gravity. These expected colors can be compared 
to the observed colors to calculate the reddening and 
extinction. The expected color is {B — V)o = —0.224 for 
a 22000 K sub-giant star, an d (B - V)o = -0.2 31 for a 
23000 K main seque nce star (jBesseh et al.lll998[) . From 
IPietsch et all (|2009[ ). the mean magnitudes are ma' = 
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Fig. 2. — Observed spectrum and best-fit model for XMMU J013236. 7+303228. The blue line is the average of the ten good spectra 
obtained on UT 2009 October 17, shifted to the rest wavelength using velocities from Table [T] The spectrum is smoothed with a 5 pixel 
(3 A) boxcar for plotting. The red line is the best-fit template spectrum with T = 22000K, log(g) = 3.5, frotsini = 250kms~^, solar 
metallicity. The template is reddened using Ay = 0.395, and scaled appropriately. Shifting the spectra to the rest frame blurs out the Ca 
II interstellar line and the 4430 A interstellar band. The observed higher Balmer lines are less strong than those of the model, suggesting 
that the surface gravity is slightly higher than logg = 3.5 (consistent with our estimated spectral type and with our fit results; see text). 

and Av = 0.425(3) with xV^OF = 1.22. The tempera- 
tures are consistent with those expected for a B1.5 sub- 
giant star, and the extinction is wit hin the range d erived 
from photometric measurements. IPietsch et al.l ()20Q9f ) 
obtained a higher extinction for the target, which ex- 
plains w hy they estima t ed th e source temperature to be 
higher. iMunari et al.l (|2005 [) templates are calculated 
in temperature steps of 1000 K in this range. For fur- 
ther analysis, we use the best-fit template: \og{g) = 3.5 
and T = 22000 K. Since log(5) is shghtly higher than 
this value, for completeness we also give results us- 
ing the best fit template for \og{g) = 4.0, which has 
T = 23000 K. For both these templates, the best-fit ex- 
tinction is Ay = 0.395(3). We then keep T and log{g) 
constant and vary v^t sini. For both the \og{g), T com- 
binations, we measure Wrot sini — 260 ± 5 km s^ . Finally, 
using the same templates but with varying metallicity, 
we get the best fits for [M/H] = 0. The 0.5 dex steps in 
[M/H] are too large to formally fit for uncertainties. 

Next, we calculate the luminosity of the object to ob- 
tain a radius - temperature relation. The distance mod- 
ulus to M33 is (w- M)M . -js = 24.54 ±0.06 (d = 809 kpc; 



21.03 ± 0.02, rrir' = 21.36 ± 0.02. Using the lJester et all 
(|200El) photometric transformations for blue, U — B < 
starfQ, mv = 21.21 ± 0.03, {B - V)ohs = -0.09 ± 0.04. 
The color excess is E{B -V) = {B- V)ohs - {B - V)o = 
0.14 ±0.04. Using the standard ratio of total-to-selective 
extinction, Ry = 3.1 we get Ay = 0.43 ± 0.12. For 
comparison, the foreground extinction to M33 is Ay = 
0.22. 

We measure various stellar parameters by fitting our 
combined, flux calibrated spectrum with model atmo- 
spheres (taking into account the instrumental broaden- 
ing; see Section 13. 2p . For a Roche lobe filling com- 
panion (discussed in Section 3]) , we expect a radius of 
6-10 i?0, surface gravity \og{g) c± 3.7 (consistent with 
our luminosity class), and projected rotational velocity 
?;i.otsini ~ 250kms~^. We use these values as start- 
ing points to sele ct mode l atmo spheres from a grid cal- 
culated by Muna rTet al.l (|2005[ ). These templates are 
calculated in steps of 0.5 dex in log(g), so we use mod- 
els with log(g) = 3.5, 4.0. For the initial fit, we as- 
sume tij-otsinz = 250 km s^^ and solar metallicity. The 
only free parameters are a normalization a nd an extinc- 
tion. We use extinction coefficients from iCoxl ()2000[ ). 
assuming Ry = 3.1. We find that the best-fit model 
for log(5) = 3.5 has T = 22100 ± 40K and Ay = 
0.401(3), with xVDOF = 1.13 for 3600 degrees of free- 
dom, while for \og{g) = 4.0, we get T = 23500 ± 50K 

6 y = 3-0.59(3-r)-0.01±0.01; B-V = 0.90(3-t-)-|-0.21±0.03. 



iMcConnachie et al.l [200a iFreedman et al.l 120011 ). which 
gives, accounting for the reddening of Ay = 0.4, My = 
—3.74 ± 0.07. The bolometric luminosity of a star 
is related to its temperature and radius by Lboi oc 
jl^rpi To obtain the visual luminosity, one must apply 
a temperatur e-dependent b olometric correction, BC = 
Mboi - My. iTorresI ((Moh give formulae for bolomet- 
ric correction as a power series in log(T). We calculate 
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Fig. 3. — Observed red spectrum for XMMU J013236. 7+303228. The blue line is the median combination of 48 spectra obtained on UT 
2009 October 17, shifted to the rest wavelength using the velocity solution from Table [2] The flux axis differs from Figure |2] The spectrum 
is smoothed with a 5 pixel (4 A) boxcar for plotting. The source is much fainter than the sky in this wavelength region. For reference, the 
green dashed line shows the extracted sky spectrum, scaled down by a factor of 100 for plotting. We do not use the red side spectra for 
any fitting. 




Fig. 4. — Radius — temperature constraint from the observed 
luminosity. We calculate the absolute visual magnitude from the 
apparent magnitude, distance modulus and best-fit extinction. We 
then apply a temperature-dependent bolometric correction to cal- 
culate the bolometric magnitude. For any given temperature in 
this range, the uncertainty in radius is about 3%. 



BC = -2.11(- 2.21) for T = 22000(23000) K (which are 
consistent with iBessell et ah] ([1998. ) tables for main se- 



quence stars.). After some basic algebra, we obtain: 



Slog 



BC{T) 



MhohQ -mv + {m- A'f )m33 + Ay 
8.48 ± 0.07 



(1) 
(2) 



Here, Mboi,© = 4.75 CBesse h et al]|1998D . The resultant 
radius-temperature relationship is shown in Figured For 
T = 22000(23000) K, we infer R = 9.1(8.7)±0.3 Rq. The 
absolute magnitude and radius are both consistent with 
a B1.5 sub-giant. 

3.2. Orbital parameters 

We measure the radi al velocities of the B star using 
model stellar spectra bv lMunari et al.l (|2005| ). as follows. 
For each observed spectrum, we measure the seeing us- 
ing the width of the spectral trace. We then generate 
an instrument response function by taking a Gaussian 
matched to the seeing, truncating it at the slit size, and 
convolving it with the pixel size. Flux calibrated syn- 
thetic spectra (templates) are convolved with this in- 
strument response, then redshifted to a test velocity. 
Then we redden the template using the measured value 
of ext inction. Ay = 0.395, with coefficients fr om iCoxl 
(|2000D . We use IDlH mpfit (jMarkwardti I2009D to cal- 
culate the reddening and normalization to match this 

^ http:/ /www. ittvis.com/ProductServices/IDL.aspx 
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spectrum with the observed spectrum, and measure the 
X^. By minimizing the over test velocities, we find 
the best-fit velocity and the error bars. This velocity 
is converted to a barycentri c radial-velocity using the 
baryvel routine in Astrolib (jLandsman' 1993). Tabled] 
lists the radial velocities for all spectra, measured using 
the best-fit stellar template: T = 22000 K, log.g = 3.5, 
[M/H] = 0.0 and Wrot sini = 250kms~\ 

Red side spectra are not useful for radial-velocity mea- 
surement for several reasons. The flux of the B star at 
redder wavelengths is lower than the blue wavelengths, 
and there are fewer spectral lines in this range. Also, 
the background noise is higher, from the large number of 
cosmic rays detected by the LRIS red side and from in- 
trinsic sky emission. Hence, all further discussion omits 
red side spectra. 

We calculate an orbital solution for the B star using 
these radial velocities. Owing to the short 1.73 d period 
of the system, we assume that the orbit must be circu- 
larized. The orbital solution is then given by: 

vit) = 7opt + Kopt sin ^27r ^ (3) 

where 7opt is the systemic velocity, i^opt is the projected 
semi-amplitude of radial-velocity, and Tq is the epoch 
of mid-eclipse. We adopt P = 1.732479 ± 000027 and 
To = 2453997.476 ± 0.006 from lPietschetall (;2009). We 
obtain 7opt = — 80±5kms^^, and iiTopt = 64±12kms~^. 
The best-fit has x^/DOF = 5.1 for 16 degrees of freedom, 
which is rather poor. We find that an additional error 
term Av — 15kms~^ needs to be added in quadrature to 
our error estimates to obtain Xred = 1- We attribute this 
to movement of the star on the slit. If the target has sys- 
tematic offset of 0".l from the slit center over the entire 
30 min exposure, it shifts the line cen troids by about 0.45 
pixels or 18kms^^. For comparison. Ivan Kerkwijk et al.l 
poll) find a similar scatter (13 kui/s~^) in their observa- 
tions of a reference star when using LRIS with a similar 
configuration (600/4000 grating, 0".7 slit). Future ob- 
servations should orient the spectrograph slit to obtain 
a reference star spectrum to correct for such an offset. 




55120.50 55121.00 55121.50 55122.00 55122.50 

HJD - 2400000 

Fig. 5. — The radial-velocity curve for 
XMMU J013236.7-I-303228. Solid blue circles denote points 
which were used in the final fit, and gray crosses are points which 
were rejected based on observing logs or a signal-to-noise ratio 
under 10 (Table [TJ- The solid blue curve is the best- fit orbital 
solution, and the dashed blue line is the mean systemic velocity. 



TABLE 2 

System parameters for XMMU ,1013236.7+303228 



Property Value 



Calculat ed using P = 1 . 73247 9±0. 000027 and To = 2453997. 476± 
0.006 from lPietsch et al.l ll2009f) . 

^ When we allow To to vary, we obtain 7opt — —80 it 5kms~^ 
and Kopt — 64 it 12kms~^. The analysis in this paper uses the 
IPietsch et al.l ll2009fl value of To . 

When we fit Equation[3]to data including the 15 km s~^ 
error in quadrature, we obtain 7opt = —80 ± 5kms~^, 
and A'opt = 63 ± 12kms"^ (Figure [3 Table [2]). For the 
epoch of observations, the uncertainty in phase is 0.019 d. 
If we allow To to vary, we get To,fit ='2453997.489±0.019 
(Heliocentric Julian Date), 7opt — — 80±5kms~^, and 
ATopt = 64± 12kms^^. These values are coii sistent with 
those obtained using the IPietsch et al.l ()2009[ ) ephemeris. 
Hence, for the rest of this paper, we simply assume their 
best-fit value for Tq. 

To investigate the sensitivity of the result to the choice 
of the stellar template, we repeat the measurement with 
a variety of templates. We vary the temperature from 
09 (33,000 K) to B3 (18,000 K) spectral classes. As be- 
fore, we use templates with log{g) ~ 3.5,4.0; Wrot sini — 
250 km s~^ and solar metallicity. Repeating the orbit cal- 
culations for each of these models, we find that the sys- 
temic velocity 7opt may change between models being fit: 
the extreme values are -77±6kms"' and -90±5kms"\ 
a 1.7(7 difference. We suspect that the reason for this 
variation is the difference in shape of the continuum, as 
the extinction was held constant in these fits. For hotter 
templates with a steeper continuum, the red side of a 
line has lower flux than the blue side, so lowest x^ will 
be obtained at a slightly higher redshift, as seen. The 
magnitude of this effect should be independent of the 
intrinsic Dopplcr shift of the spectrum, and should not 
affect the velocity semi-amplitude ATopt. This is indeed 
the case: A'opt is constant irrespective of templates. The 
extreme values are 62±12kms^^ and 63±12kms~^, dif- 
fering by less than O.lu. Dynamical calculations depend 
only on ATopt, hence they are robust to the selection of 
template. 

4. COMPONENT MASSES 

The general method for accurately determining masses 
in X-ray binaries (jJoss fc Rappaport..l984) requires mea- 



Prom Pietsch et al. (2009) 
Period (P) 

HJD of mid-eclipse (To) 
Eclipse half-angle (de) 

This work 

Systemic velocity'' (7opt) 
Velocity semi-amplitude^ (fcopt) 
HJD of mid-eclipse'" (To) 

Spectroscopically inferred 

OB star spectroscopic mass (Afopt) 
NS mass (Afx) 

Distance-based calculations 
OB star mass (Mopt) 
NS mass (Mx) 



1.732479 
2453997.476 
30°. 6 



63 

2453997.489 



11 
2.0 



18.1 
2.7 



±0.000027 

±0.006 

±1°.2 



±5kms 

±12kms-i 

±0.019 



±1M0 
±0.4 M0 
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suring the orbit for both components, as well as having a 
constraint on the orbit inclination. In general, the mass 
(Ml) of a component is expressed in terms of five pa- 
rameters: the orbital period P, the radial-velocity semi- 
amplitude of the companion {K2), the eccentricity e, the 
orbital inclination i and mass ratio q = M1/M2. The 
first three parameters can be readily obtained by char- 
acterizing the orbit of either component through pulse 
timing of the NS in the X-ray or radio, or by spectroscop- 
ically measuring the radial-velocity of the optical com- 
panion at optical or infra-red wavelengths. Determining 
the mass ratio requires measuring orbital parameters for 
both components: q = K2/K1. In an eclipsing system, 
the inclination can be constrained to be nearly edge-on 
{i « 90°), with a lower limit derived from eclipse dura- 
tion and Roche-lobe arguments. Using these measure- 
ments, masses of both components in the system can be 
directly determined without mode l assumptions (see for 
example, Ivan der Meer et"aLll2007HMason et al.ll2Qll[ ). 

Because no pulsations have been detected from 
the compact obje ct in XMMU J013236. 7+303228 
(jPietsch et all 120091 ) . we need one more constraint in 
addition to the radial-velocity semi-amplitude of the 
donor. In Section 14.11 we use the spectroscopically 
inferred mass of the primary to calculate mass of the 
secondary from the mass ratio. Because we know the dis- 
tance to M33, and hence to XMMU J013236. 7+303228, 
we can also estimate the physical size of the secondary 
from the distance, its luminosity and temperature. 
This provides a cross check on the mass determination 
derived from the spectroscopic donor mass (Section l4.2p . 

4.1. The spectroscopic method 

In the following, we will denote the masses of the com- 
pact object and donor star by Mx and Mopt respectively. 
The mass of the compact object (Mx) is related to the 
radial-velocity of the B star (i^opt) as follows: 



Mx = 



^gpt^(l-e ^) 
2nG sin^ 



2\3/2 



1 



(4) 



where q = Mx/A/opt is the ratio of masses, defined so 
that higher values of Afx relate to higher values of q. P 
is the orbital period of the binary, e is the eccentricity, 
and i is the inclination of the orbit. For eclipsing systems, 
the inclination is constrained by. 



sm^ = 



cos Or, 



(5) 



where i?L is the volume radius of the Roche lobe, 
a is the semi-maj or axis, and /3 i s th e Roche 
lobe fining factor poss fc Rappapori] I1984D . For 
XMMU JO 132 36. 7+303228 the eclipse half-angle is 9^ = 
30°.6± 1°.2 ( Pietsch et al.ll2009D . Owing to the short or- 
bital period, we assume that the orbit is circular and the 
B star rotation is co mpletely synchron ized with its orbit. 
For co-rotating stars. lEggletonl (|19830 expresses Ri,/a in 
terms of q: 



0.49g 



-2/3 



0.6g-2/3 + ln(l + g-l/3) 



(6) 



The constant, relatively high X-ray luminosity, sus- 
tained over the non-eclipsed parts of the orbit, strongly 



indicates that accretion is occurring via Roche lobe over- 
flow. In Roche lobe overflow, matter flowing through the 
Lagrangian point may form a disc around the compact 
object before being accreted onto it. This disc can occult 
the compact object, causing periods of low X-ray lumi- 
nosity. Both these characteristics are seen in the X-ray 
lightc urves of XMMU J013236.7+303228 ()Pietsch et al.l 
120091 ). If mass is being accreted onto an object in a spher- 
ically symmetric manner, the accretion rate is limited by 
the Eddington rate, Mgcid and the peak luminosity is 
^Edd/i© = 3 X 10^M/Mq. For a 1.5-2.5 Mg compact 
object, LEdd = 1.8 - 3 X lO^^ergs"^ At its bright- 
est, the source luminosity in the 0.2-4.5keV band was 
2.0 X lO^^ergs"^ — about O.lLEdd- This luminosity was 
sustained throughout Chandra ObsID 6387, whic h cov- 
ered about 0.72 d of the non-eclipsed orbit (Piet sch et al.l 
120091 ). Comparable flux was observed in the non-eclipsed 
parts of the orbit (0.73 d) in Chandra ObsID 6385. Such 
high luminosity sustained over significant parts of the or- 
bit is not observed in wind-fed systems, which have typi- 
cal luminosities an order of magnitude smaller. Further, 
the short 1.73 d orbital period is not consistent with Be 
X-ray binary or wind- fed systems. We conclude therefore 
the B star fills its Roche lobe. 

Radius / Mass for the OB star 
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Mass (Mq) 

Fig. 6. — Mass-Radius relation for the OB star. The sohd hnes 
show the Roche lobe radius for the OB star, assuming a 2.4 Mq 
(upper blue line) and a 1.4 Mq (lower red line) neutron star. The 
dashe d green lines are con tours for log(g) in steps of 0.5 dex, for 
which iMunari et al.l II2005I ) synthetic spectra are available. In the 
range 3.5 < log((7) < 4.0, dotted \og{g) contours are separated by 
steps of 0.1 dex. The stars denote the masses and radii of typical 
isolated B stars. 



The Roche lobe radius as a function of B star mass is 
plotted in Figure [S] for various neutron star masses. For 
a Roche-filling companion, we infer 3.6 < \og{g) < 3.8. 
For a NS mass in the range 1.4 - 2.4 M© and B star 
mass 8 - 20 M©, the B star radius lies in the range 6 - 
lOi?0. Further, the assumed synchronous rotation re- 
quires the surface rotational velocity to be in the range 
200kms~^ < Wrot ^ 285kms~^. These values are con- 
sistent with those derived in Section [3+1 

Figure [7] plots Mx as a function of Afopt, calculated by 
solving Equations min] and [5] under known constraints. If 
we know Mopt, we can calculate Mx. We use the physical 
properties of the primary (Section 13. 1[) to estimate the 
mass of the primary by comparing it with stellar evolu- 
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4 8 12 16 20 

Fig. 7. — The solid blue line shows the compact object mass (Afx) as a function of OB star mass (Mopt). Dashed blue lines are itl-o" 
errors. The estimated mass of the B1.5IV primary is 11 Mq, marked by the vertical dotted line. The corresponding mass of the neutron 
star is 2.0 ± 0.4 M©. 




Fig. 8. — Evolutionary tracks and isochrones for massive stars (adapted from lBrott et ani2011l ). Left panel: the conventional HR diagram 
with luminosity and temperature. The primary star has 22000 K < T < 2300 K. We calculate luminosity from the observed mv, and 
allow a 0.5 mag offset to calculate the lower limit (see discussion in Section l4.1l l. This region is shown by a shade d gr ay box. Right panel: 
same, but plotted as \og{g) versus T. log((/) for the primary is constrained from Roche lobe arguments fSection I3.1[l . From both panels, 
we see that the primary is consistent with a 11 Mq, 15Myr object. 
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tionary models, assuming that binary evolution has not 
drastically changed the mass-luminosity relation. First, 
we place t he primary o n a HR diagram (Figure ^ using 
models by iBrott et al.l ()2011| ) . We conservatively allow 
for a 0.5 mag error in luminosity. We also plot evolution- 
ary tracks on a log{g) - T figure, to utilize the stricter 
constraints on \og{g) from Roche lobe arguments. From 
these plots, we see that the primary is approximately a 
15Myr old, ~ 11 Mq star. T his is cons istent with the 
typical mass of a B1.5IV star (|Co>dl2000( ). 

For Mopt = 11 Mq, we calculate Mx = 2.0 ± 0.4 M©. 
From evolutionary tracks, we estimate that the uncer- 
tainty in Mopt is ~ IMq (Figure [S]), corresponding to 
AMx = 0.12, much smaller than the uncertainty aris- 
ing from Ai^opt. Adding this in quadrature with the 
uncertainty in the Mx - Mopt conversion, we conclude 
Mx = 2.0±0.4Mo. 

4.2. Masses from Roche lobe constraints 

In Section 13.11 we calculated the radius of the primary 
from its apparent magnitude, temperature and the dis- 
tance to M33 (Equation [Ij . Since the primary is filling 
its Roche lobe, the stellar radius is equal to the Roche 
lobe radius (Rh). This additional constraint can be used 
in Equations H]- [6] to solve for Mx and Mopt- 

We calculate the probability density function (PDF) of 
component masses as follows. For every pair of assumed 
masses (Mx, Mopt), we use the period P to calculate 
the semi-major axis a. Then we calculate Ri^/a from 
Equation [6] and substitute it in Equation [5] to calculate 
sini. Using P, a and sini we calculate the expected 
semi-amplitude of the radial-velocity: 

27rasinz Mx , . 

' " P (Mx + Mopt) ^' 

Next, we calculate the probability for obtaining a cer- 
tain value of i?L and K2, given the measured radius R 
(Section [3TT|) and Kopt (Section [3^ : 

P{RuK,) ^ exp ) exp (^-^-^^ j 

(8) 

Here, we are making a simplifying assumption that the 
Roche volume radius (Equation|6l) is same as the effective 
radius from photometry (Equation [1]) . We convert this 
PDF to a probability density as a function of Mx, Mopt 
by multiplying by the Jacobian 9(i?, K^pt) / d{Mx, Mopt). 

The results are shown in Figure [S] Red contours show 
the 68.3% and 95.4% confidence intervals for masses for 
the best-fit template (log(.g) = 3.5 and T = 22000 K). 
The panel on the left shows the PDF for Mx marginal- 
ized over Mopt. Similarly, the lower panel shows the 
PDF for Mopt, marginalized over Mx- In these pan- 
els, the the solid, dashed and dotted lines show the peak 
and 68.3%, 95.4% confidence intervals respectively. We 
obtain Mx = 2.7lJ^;^M0, and Mopt = IS.ltf^MQ- 
For completeness, fits for the less likely scenario with 
\og{g) = 4.0 and T = 23000 K are shown in blue. In this 
case, Mx = 2.5 ± 0.6 Mg, Mopt = 16.ltJ:?MQ. 

To test the validity of this technique, we apply it to 
two well-studied targets: LMCX-4 and SMCX-1. We 
find that the basic application of our method is over- 
estimating the mass (Table[3]). Part of this discrepancy is 



likely related to equating the Roche volume radius to the 
effective photometric radius (Equation [8]). The Roche 
volume radius (Equation |6]) is the radius of a sphere with 
the same volume as the Roche lobe of the star. The 
effective photometric radius (Equation [T]) is the radius of 
a sphere with the same surface area as the star. Since 
an ellipsoid has a larger surface area than a sphere of 
the same volume, the actual Roche volume radius will 
be smaller than the effective photometric radius. Using 
a larger radius for the Roche lobe increases the masses 
of both the components in the binary. 

Putting it another way, if we use the Roche volume 
radius to calculate the brightness of the star, we will get 
a number lower than the o bserved brightness. A similar 
discrepancy is observed by iMassev et al.l (|2012[ ) in mas- 
sive binaries in the LMC. They find that the absolute 
magnitude of LMC 172231 calculated using a spherical 
approximation is 0.45 mag fainter than observed, while 
for the triple system [ST92]2-28, the numbers are con- 
sis tent within errors. Using s ystem parameters derived 
by Ivan der Meer et all ()2007D . we find a similar offset 
of 0.45 mag for LMCX-4 and 0.2 mag for SMCX-1. If 
we incorporate this uncertainty by allowing offsets of - 
0.4 mag, we get Mx = 2.2^° J Mq and Mopt = 13±4M0. 

Thus while this method has potential, more detailed 
modeling of the primary is clearly required to accurately 
infer component masses. 

5. CONCLUSION 

From our spectroscopic measurements we find that the 
donor star in XMMU J013236.7-K303228 is a B1.5IV sub- 
giant with effective temperature T = 22000 - 23000 K. 
Assuming a circular orbit, we measure a mean sys- 
temic velocity 7opt = —80 ± 5kms~^ and velocity semi- 
amplitude i^opt = 63 ± 12kms~^ for the B star. M33 
is nearly face-on, with reces sion velocity of —179 km s^^ 
(|de Vaucouleurs et al.lll99l"l ) - so this binary seems to be 
moving away from the disc at lOOkms"^. 

Using the physical properties of the B star determined 
from our optical spectroscopy we find a mass for the 
donor of Mopt = 11 ± 1 Mq. This mass is based on stel- 
lar evolution models, and will be reasonably accurate so 
long as binary evolution has not significantly altered the 
mass-luminosity relation. However, it is difficult to test 
this assumption based on any available observations. Us- 
ing this spectroscopic mass, we calculate the mass of the 
compact object, Mx — 2.O±O.4M0. This is higher than 
the canonical 1.4 Mq for neutron stars, but comparable 
to masses of other neutron stars in X-ray binaries such 
as the HM XB Vela X-1 (1 .88 ± 0.13 Mq; iBarziv et all 
2001; Quai ntrell et al.l 120031 ) or the Low Mass X-ra y bi- 
naries Cyg X-2 (1.71±0.21 M r7,:lCasares et al."2010') and 
4U 1822-371 (1.96±0.35MQ: lMunozDarias et al...2005) . 
Since no pulsations have been detected we have only in- 
direct evidence, based on the X-ray spectrum, that the 
compact object is a neutron star. However, the mass we 
derive here is smaller than would be expected for a black 
hole. 

Based on the stable X-ray flux, we infer that the donor 
is transferring mass to the neutron star by Roche lobe 
overflow. By equating the Roche lobe radius to physical 
radius of i? = 9.1 ± 0.3 i?Q, derived from the known dis- 
tance to M33, combined with the stellar luminosity and 



10 



Bhalerao et al. 




- 












_ 




/V / 




y ; 


- 


U~)l 
CT)/ 
/ 


u / 
oV / 


/ / h 


/ ^ 
















1 












0.6 0.8 5 

0.25 F 



0.05 : 
0.00 E 




Fig. 9. — Probability density plot for neutron star mass (Mx) as a function of OB star mass (Mopt)- Red contours show the 68.3% 
and 95.4% confidence intervals for masses using the best-fit template (log(g) = 3.5 and T = 22000 K). The panel on the left shows the 
PDF for Afx marginalized over Mopt. Similarly, the lower panel shows the PDF for Mopt, marginalize over Mx. In these panels, the 
the solid, dashed and dotted lines show the peak and 68.3%, 95.4% confidence intervals respectively. We obtain Mx = 2.7 



+0.7 



0.6^^0. 
In this case. 



Mopt = 18.1_]^'gM0. Contours for the the less likely scenario with log(g) = 4.0 and T = 23000 K are shown in blue. 
Mx = 2.5±O.6M0, Mopt = 16.lj:i;7MQ. 



temperature, we derive an additional orbital constraint. 
"From a first pass calculation with a spherical approx- 
imation for the shape of the primary, we find substan- 
tially larger masses. However, applying this technique 
to the well-studied binaries LMC X-4 and SMC X-1, 
both of which have measured component masses, we find 
it consistently overestimates the compact object mass. 
This is likely because the Roche surface is not spher- 
ical but elongated, which is not taken into account in 
our estimate. Attempting to account for that, we infer 



Mx = 2.7 



+0.7 
-0.6 



Mq and Mopt = 18.1 ± 2.0 Mq. Future ef- 
forts to more accurately model the system geometry will 
improve the accuracy of this technique, which is applica- 
ble generally to Roche lobe overflow systems with known 
distances. 
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TABLE 3 

System parameters for SMCX-1 and LMCX-4 



Property 
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LMCX-4 


Period (P) 
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1.41 d 
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zb — oO.o 


2 1 ±2 
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20.2±l.lkms-i 


35.1±1.5kms-i 


Companion spectral type 


BOi 


08111 


Companion T^ff 


ononri 

zyouo K 


OOOOO K 


Distance (d) 


eo.eiikpc'^ 


49.4 ±1 kpc'" 


Visual magnitude {my)'^ 
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Extinction (Ay)'^ 
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NS mass (Mx) 
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Our calculation 
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OB star mass (Mopt) 






van der Meer et al. f2007) 


« 15. 7 Me 
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Our calculation 


21.1 1^.2 M0 


29.1 ttlMQ 


Note. — Data from Ivan der Meer et al.l (120071). 
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IFreedman et al.l (I200l]l. 
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Foreground extinction to galaxy only. Error bars arc approximate, tlic dominant 
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